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ABSTRACT 


We utilize recent NuSTAR observations (co-added depth ~ 55-120 ks) of PG 1001 +054, PG 1254+047, and 
PHL 1811 to constrain their hard X-ray (Z; 5 keV) weakness and spectral shapes, and thus to investigate the 
nature of their extreme X-ray weakness. These quasars showed very weak soft X-ray emission, and they were 
proposed to be intrinsically X-ray weak, with the X-ray coronae producing weak continuum emission relative 
to their optical/UV emission. However, the new observations suggest an alternative explanation. The NuSTAR 
3—24 keV spectral shapes for PG 1001 +054 and PHL 1811 are likely flat (effective power-law photon indices 
Tere = Low and Deg = pas, respectively), while the shape is nominal for PG 1254 - 047 (Veg = 1.8 3: 0.3). 
PG 1001+054 and PHL 1811 are significantly weak at hard X-ray energies (by factors of ~ 26—74 at rest-frame 
8 keV) compared to the expectations from their optical/UV emission, while PG 1254 4-047 is only hard X-ray 
weak by a factor of ~ 3. We suggest that X-ray obscuration is present in all three quasars. We propose that, as 
an alternative to the intrinsic X-ray weakness + X-ray obscuration scenario, the soft and hard X-ray weakness 
of these quasars can be uniformly explained under an obscuration-only scenario. This model provides adequate 
descriptions of the multi-epoch soft and hard X-ray data of these quasars, with variable column density and 
leaked fraction of the partial-covering absorber. We suggest that the absorber is the clumpy dust-free wind 
launched from the accretion disk. These quasars probably have super-Eddington accretion rates that drive 
powerful and high-density winds. 
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1. INTRODUCTION the X-ray and UV emission from large samples of radio- 
quiet AGNs have revealed that the X-ray flux is closely cor- 
related with the optical/UV flux, indicating a strong phys- 
ical connection between the accretion disk and the X-ray 
corona. The correlation is typically expressed as a nega- 
tive relation between the X-ray-to-optical power-law slope 


Active galactic nuclei (AGNs) generally produce luminous 
X-ray emission, which is believed to originate largely from 
the accretion-disk corona in the vicinity of the central su- 
permassive black hole (SMBH) via inverse-Compton scatter- 
ing of the optical/UV seed photons from the accretion disk 


(e.g., Turner & Miller 2009; Done 2010; Gilfanov & Merloni parameter (aox)' and the 2500 A monochromatic luminos- 
2014; Fabian et al. 2017). AGN X-ray continua typically VE goo 20d it is highly significant across a broad pop- 
have an intrinsic power-law shape (Ng c ET), and the mean 

value of the photon indices (I) for radio-quiet AGNs is lagy is defined as aox = 03838108 firey Iso A) (Tananbaum et al. 


around 1.9-2.0 with a scatter of ~ 0.2 (e.g., Reeves et al. 1979), where fxev and f,- & are the rest-frame 2 keV and 2500 A flux 
1997; Just et al. 2007; Scott et al. 2011). Observations of densities. 
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ulation of AGNs, ranging from moderate-luminosity AGNs 
to the most luminous quasars (e.g., Strateva et al. 2005; 
Steffen et al. 2006; Just et al. 2007; Lusso & Risaliti 2017; 
Liu et al. 2021). 

The observed X-ray emission from AGNs may be modi- 
fied by line-of-sight obscuration, resulting in lower observed 
X-ray fluxes than those expected from the aox—L * re- 
lation. A common approach to parameterize the amount of 
X-ray weakness uses the Aaox parameter, defined as the dif- 
ference between the observed oox value and the aox value 
expected from the aox—L A relation; Aaox = —0.3838 
thus corresponds to a factor of X-ray weakness of 10 at 
rest-frame 2 keV. Type 2 AGNs are generally X-ray ob- 
scured, likely due to the dusty “torus” (e.g., Netzer 2015; 
Hickox & Alexander 2018). Type 1 AGNs may also have 
X-ray obscuration from largely dust-free gas.^ For exam- 
ple, broad absorption line (BAL) quasars, which are char- 
acterized by blueshifted broad UV absorption lines (e.g., 
C Iv A1549), generally show weak X-ray emission (e.g., 
Gallagher et al. 2002, 2006; Fan et al. 2009; Gibson et al. 
2009). One frequently adopted physical model for BAL 
quasars is the disk-wind model, where the observed BALs 
originate from an outflowing equatorial wind launched from 
the accretion disk and radiatively driven by UV-line pressure 
(e.g., Murray et al. 1995; Proga et al. 2000; Matthews et al. 
2016). This model usually invokes “shielding” gas between 
the wind and nucleus or a clumpy wind (e.g., Baskin et al. 
2014; Matthews et al. 2016; Giustini & Proga 2019) to pro- 
vide obscuration of the nuclear extreme UV (EUV) and 
X-ray radiation which might otherwise overionize the wind 
and hamper radiative acceleration. BAL quasars are con- 
sidered generally to have larger inclination angles than non- 
BAL quasars, with the line of sight to the UV continuum 
region of the accretion disk intersecting the wind, leading 
to the observed BALs. The line of sight to the X-ray emit- 
ting corona, though not necessarily the same as the UV 
line of sight, is likely also through the shielding gas or the 
clumpy wind, resulting in the often observed X-ray weakness 
(e.g., Figure 1 of Luo et al. 2013). Besides BAL quasars, 
a small fraction (5.8 + 0.7%) of non-BAL type 1 quasars 
have been found to be X-ray weak, likely due to absorp- 
tion (e.g., Pu et al. 2020). They may share a similar nature 
to the BAL quasars; they do not show any UV BALs prob- 
ably due to geometric effects (e.g., small inclination angles) 
or a low velocity of the wind along the UV line of sight (e.g., 
Giustini & Proga 2019). 

A small subset of AGNs has been proposed to be intrin- 
sically X-ray weak, producing much less X-ray emission 


? Similar obscuration from dust-free gas might also be present in some of the 
type 2 AGNs, though usually not distinguishable from the torus obscura- 
tion. 


than expected from the aox—L o relation. These candi- 
dates are observed to be significantly X-ray weak with no 
clear evidence of X-ray obscuration. X-ray weakness caused 
by X-ray obscuration is usually identified from an X-ray 
spectral shape that is flatter than the intrinsic I' ~ 2 power 
law, as soft X-ray photons are more heavily absorbed than 
hard X-ray photons. The effective power-law photon in- 
dex (Tere) of an obscured spectrum should be smaller than 
2; in case of heavy absorption, the 0.5-5 keV Teg can even 
reach a negative value. However, due to the frequent ap- 
pearance of partial-covering absorption in AGN X-ray spec- 
tra (e.g., Immler et al. 2003; Ricci et al. 2017; Leighly et al. 
2019) with a small fraction of the intrinsic coronal emission 
leaking through the absorber, a soft X-ray (<5 keV) spec- 
trum alone might be insufficient for identifying heavy (Nu = 
5 x 10? cm?) or Compton-thick (Ny > 1.5 x 10% cm?) 
X-ray obscuration, as the X-ray emission could be extremely 
weak (e.g., Aaox < —0.3838 or X-ray weakness factor of 
> 10) yet the spectral shape is nominal with Ief z 2, dom- 
inated by the leaked component (e.g., see Section 4.2 below 
for illustration). On the other hand, the hard X-ray (> 5 keV) 
spectrum in this case should generally still be flat with Date © 
0-1, as the Compton-reflection component, from either the 
absorber or other reflectors (e.g., disk or torus), is expected to 
dominate (e.g., George & Fabian 1991; Comastri et al. 2011; 
Gandhi et al. 2014; Rovilos et al. 2014). It is based on these 
arguments that a small number of quasars have been pro- 
posed to be intrinsically X-ray weak; they are significantly 
X-ray weak yet with nominal (Ief œ~ 2) hard X-ray spec- 
tral shapes. The low-redshift (z < 1) candidates include a 
few BAL quasars that have NuSTAR observations (Luo et al. 
2013, 2014), and the high-redshift (z ~ 1.5-3) candidates 
include a few Large Bright Quasar Survey (LBQS) BAL 
quasars with Chandra observations (Liu et al. 2018) and 
a few luminous quasars with XMM-Newton observations 
(Nardini et al. 2019). Given the significant X-ray weakness, 
the X-ray spectra of these candidates have very limited pho- 
ton statistics, and thus the spectral shapes are poorly con- 
strained and sometimes even rely upon stacking analysis to 
obtain information on their average spectral properties (e.g., 
Luo et al. 2014). 

One quasar that is considered a prototypical example of 
intrinsically X-ray weak quasars is PHL 1811, a non-BAL 
quasar at z = 0.192 (Leighly et al. 2007b). Itis X-ray weak by 
a factor of ~ 40-130 with a steep 0.3-5 keV XMM-Newton 
spectrum (I =~ 2.3). There are no published hard X-ray 
constraints, but its fast X-ray variability (flux varying by 
a factor of ~ 4 in 12 days) argues against a heavily ob- 
scured spectrum dominated by distant reflection/scattering 
(Leighly et al. 2007a,b). PHL 1811 has distinctive UV 
emission-line properties, including a small C Iv rest-frame 
equivalent width (REW), a large C Iv blueshift, and strong 


UV Fe II and Fe III emission. Luo et al. (2015) conducted 
a systematic survey of the X-ray properties of PHL 1811 
analogs, a sample of high-redshift quasars selected with UV 
emission-line properties similar to those of PHL 1811. This 
sample of 18 PHL 1811 analogs turned out to be almost ex- 
clusively (17/18; 94%) X-ray weak, by an average X-ray 
weakness factor of 39. However, unlike PHL 1811 itself, 
many of these PHL 1811 analogs appear to be X-ray ob- 
scured, as the sample stacked effective photon index (Teff = 
1.16732) indicates a flat spectral shape on average. Luo et al. 
(2015) speculated that Occam's razor would favor a uniform 
explanation for the X-ray weakness of all these objects and 
that perhaps hard X-ray data of PHL 1811 would reveal a 
highly obscured component. 

In this paper, we present improved hard X-ray constraints 
from deeper NuSTAR observations of two low-redshift in- 
trinsically X-ray weak quasar candidates in Luo et al. (2014), 
PG 1001+054 (hereafter PG 1001) and PG 1254+047 (here- 
after PG 1254). We also provide for the first time hard 
X-ray constraints for PHL 1811 using archival NuSTAR and 
XMM-Newton observations. Based on the results, we sug- 
gest that X-ray obscuration is present in all three quasars, 
and we propose that intrinsic X-ray weakness is not required 
to explain the X-ray weakness of these objects. The NuS- 
TAR, Chandra, and XMM-Newton data of these quasars can 
be uniformly explained with an obscuration scenario where 
the partial-covering absorber has variable column density and 
leaked fraction (partial-covering fraction). The paper is orga- 
nized as follows. We describe the soft and hard X-ray obser- 
vations and data analyses in Section 2. We present X-ray 
and multiwavelength properties of the three quasars in Sec- 
tion 3. In Section 4, we propose that these quasars are likely 
X-ray obscured and describe how an obscuration scenario 
may explain their soft and hard X-ray weakness, without in- 
voking intrinsic X-ray weakness. The absorber is likely the 
clumpy dust-free wind launched from the accretion disk. We 
summarize in Section 5. In the Appendix, we present new 
hard X-ray constraints from a recent XMM-Newton observa- 
tion of LBQS 1442—0011, a high-redshift intrinsically X-ray 
weak quasar candidate in Liu et al. (2018); the results are 
consistent with our proposed obscuration scenario. 

Throughout this paper, we use a cosmology with 
Ho = 67.4 km s! Mpc!, Qm = 0.315, and Qa = 0.685 
(Planck Collaboration et al. 2020). Uncertainties are quoted 
at a lo confidence level, and limits are at a 90% con- 
fidence level. The energy ranges used in our photomet- 
ric and spectroscopic analyses are 0.3-8 keV for Chan- 
dra observations, 0.3-10 keV for XMM-Newton observa- 
tions, and 3-24 keV for NuSTAR observations, unless oth- 
erwise specified. Due to the X-ray weakness of the three 
quasars, all X-ray spectra were grouped with at least one 
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count per bin, and the W statistic? of XSPEC (v12.10.1; 
Arnaud 1996) was used in spectral fitting. Galactic absorp- 
tion (HIAPI Collaboration et al. 2016) was included in the 
X-ray spectral modeling. 


2. OBJECT PROPERTIES, X-RAY OBSERVATIONS, 
AND X-RAY DATA ANALYSES 


2.1. Basic Object Properties 


The basic properties of the three quasars are listed in 
Table 1. PG 1001 is at z 2 0.161 with a B-band magni- 
tude of 16.1. Its HG-based single-epoch virial SMBH mass 
(Mg) is ~ 5.5 x 10’Mo, and the estimated Eddington ra- 
tio is = 0.5 (Shen et al. 2011). The full width at half max- 
imum (FWHM) of its H8 emission line is 1740 km s^! and 
thus it was classified as a narrow-line type 1 quasar (NLQ1; 
Wills et al. 2000), which refers to type 1 quasars with narrow 
(FWHM < 2000 km s^!) Hf emission lines (e.g., Footnote 
1 of Wills et al. 2000). It also shows weaker [O nr] A5007 
emission compared to typical quasars (REW smaller by a fac- 
tor of ~ 2; Vanden Berk et al. 2001; Shen et al. 2011). 

PG 1254 is a luminous quasar at z = 1.026 with a B- 
band magnitude of 15.8. The Mg I1-based single-epoch virial 
SMBH mass is ~ 4.8 x 10°Mo, and the estimated Eddington 
ratio is © 0.3 (Shen et al. 2011). There are no H£ or [O m] 
measurements available in the literature. 

PHL 1811, at z 2 0.192, is very bright with a B-band mag- 
nitude of 13.9. Its HS-based single-epoch virial SMBH mass 
is © 1.8 x 108Mo, and its estimated Eddington ratio is ~ 1.6 
(Leighly et al. 2007b). We caution that the virial SMBH 
mass and the estimated Eddington ratio are subject to large 
uncertainties, especially in the super-Eddington regime, as 
the virialization assumption may no longer be valid when 
the broad emission-line region (BELR) is likely exposed to 
large and anisotropic radiation pressure (e.g., Marconi et al. 
2008, 2009; Netzer & Marziani 2010). PHL 1811 was clas- 
sified as a NLQI with no apparent [O 11] 45007 emission 
(Leighly et al. 2007a). 


2.2. NuSTAR Observations and Data Analysis 


PG 1001 and PG 1254 were observed by NuSTAR in 2013 
with cleaned exposure times of 19.6 ks and 29.4 ks, respec- 
tively. They were not detected in the 8—24 keV band indi- 
vidually, but they were considered good candidates for in- 
trinsically X-ray weak quasars from X-ray stacking analysis 
(Luo et al. 2014). We obtained much deeper NuSTAR obser- 
vations of PG 1254 in 2019 and PG 1001 in 2020 with expo- 
sure times of ~ 100 ks (PI: W. N. Brandt), with the aim of 
detecting them individually in the 8—24 keV band and provid- 
ing improved spectral-shape constraints. PHL 1811 has an 


3 https://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/XSappendixStatistics.html. 


Table 1. Basic Object Properties and List of Observations 


L / Lgaád log L 


Object z mg  logMsnu A 
Name (Mo) (erg s7! Hz!) 
(1) Q) (3) (4) (5) (6) 
PG 1001+054 0.161 16.1 7.74 0.5 29.9 
PG 1254+047 1.026 15.8 9.68 0.3 31.5 
PHL 1811 0.192 13.9 8.25 1.6 30.9 


Nu.Gal Observatory Observation Obs. Date Exp 
(10? cm?) ID (ks) 
(7) (8) (9) (10) (11) 
2.38 XMM-Newton 0150610101 2003 May 4 8.6 
Chandra 11852 2010 Jan 11 4.9 
NuSTAR 60001122002 2013 Jun 28 19.6 

NuSTAR 60501014001 2020 May 23 101.1 
2.03 Chandra 832 2000 May 29 36.0 
NuSTAR 60001123002 2013 Jun27 29.4 
NuSTAR 60401013002 2019 Jun8 . 88.0 
4.22 Chandra 2975 2001 Dec 5 9.3 
Chandra 2985 2001 Dec 17 9.8 
XMM-Newton 0204310101 2004 Nov 1 23.5 
XMM-Newton 0761910201 2015 Nov29 39.9 
NuSTAR 60101004002 2015 Nov28 54.7 


NOTE— Cols. (1) and (2): object name and redshift. Col. (3): B-band magnitude. Cols. (4) and (5): single-epoch virial BH mass and 
Eddington ratio from Shen et al. 2011 (for PG 1001 and PG 1254) or Leighly et al. 2007b (for PHL 1811). Col. (6): 2500 À monochromatic 
luminosity from Shen et al. 2011 (for PG 1001 and PG 1254) or Leighly et al. 2007b (for PHL 1811). Col. (7): Galactic neutral hydrogen 
column density (HI4PI Collaboration et al. 2016). Cols. (8) and (9): observatory and observation ID. Col. (10): observation start date. Col. 
(11): cleaned exposure time; for NuSTAR observations, it is the average value of the FPMA and FPMB exposure times. 


archival NuSTAR observation in 2015 with a cleaned expo- 
sure time of 54.7 ks, simultaneous to a 39.9 ks XMM-Newton 
observation (PI: K. Leighly). The details of the NuSTAR ob- 
servations are listed in Table 1. 

For NuSTAR data reduction, we used HEASoft (v6.29) 
and the NuSTAR Data Analysis Software (NuSTARDAS; 
v2.1.1) with NuSTAR CALDB 20210728. We used the 
NUPIPELINE script to generate calibrated clean event files 
from the unfiltered event files of the two focal plane mod- 
ule detectors (FPMA and FPMB). For the 2019 NuSTAR ob- 
servation of PG 1254, background event rates were slightly 
elevated around the South Atlantic Anomaly (SAA), and we 
followed the recommendations from the NuSTAR instrument 
team and applied an additional SAA filter (SAAMODE = OP- 
TIMIZED, TENTABLE = YES, and SAACALC = 1) during the 
NUPIPELINE processing. We created NuSTAR images us- 
ing the Chandra Interactive Analysis of Observation (CIAO; 
v4.13)* tool DMCOPY in three energy bands: 3-24 keV (full 
band), 3-8 keV (soft band), and 8—24 keV (hard band). 

For each NuSTAR observation, we co-added the FPMA 
and FPMB images to improve sensitivity, which helps in 
detecting and characterizing faint NuSTAR sources (e.g., 
Lansbury et al. 2017). The co-added images were created 


^ http://cxc.harvard.edu/ciao/. 


by combining the FPMA and FPMB images in each of the 
three bands using the HEAsoft tool XIMAGE. These im- 
ages were then used for source detection and aperture pho- 
tometry. For each co-added image, we searched for X-ray 
sources using the CIAO tool WAVDETECT (Freeman et al. 
2002) with a false-positive probability threshold of 10? and 
wavelet scales of 2, 2.83, 4, 5.66, 8, 11.31, and 16 pixels 
(e.g., Luo et al. 2014); the NuSTAR pixel size is 2.46". In 
the 2013 observations, PG 1001 and PG 1254 were detected 
in the soft and full bands but not in the hard band. They were 
detected in all three bands in the latest deeper observations. 
PHL 1811 was also detected in all three bands. In the follow- 
ing analysis, we adopted the full-band WAVDETECT positions 
as the X-ray positions. The offsets between the X-ray and op- 
tical positions range from 3.4" to 9.4", with are typical for 
faint NuSTAR sources (e.g., Lansbury et al. 2017). 

We extracted source and background spectra using the 
NUPIPELINE script. We used 35"-radius circular source re- 
gions centered on the X-ray positions, and annular back- 
ground regions centered on the X-ray positions with inner 
radii of 120" and outer radii of 180". We verified that there 
are no sources in the background regions. For each observa- 
tion, we merged the FPMA and FPMB source spectra, back- 
ground spectra, and response files using the HEASoft tool 
ADDSPEC. 


2.3. Archival Chandra and XMM-Newton Observations and 
Data Analyses 


The archival Chandra and XMM-Newton observations of 
the three quasars are listed in Table 1? For Chandra data 
reduction, we used the CIAO CHANDRA_REPRO script to 
generate new level 2 event files. Background flares were 
then filtered with the DEFLARE script using an iterative 3c 
clipping algorithm. We created 0.5-8 keV images from the 
cleaned event files using DMCOPY. We searched for X-ray 
sources in the X-ray images using WAVDETECT with a false- 
positive probability threshold of 10% and wavelet scales of 
1, 1.414, 2, 2.828, 4, 5.656, and 8 pixels. The three quasars 
were all detected within 0.22"—0.42" of the optical positions. 
For each observation, the source spectrum was extracted us- 
ing the SPECEXTRACT tool from a circular region centered 
on the X-ray position with a radius of 2". The background 
spectrum was extracted from an annular region centered on 
the X-ray position with an inner radius of 6" and an outer 
radius of 10"; we verified that the background regions do not 
contain any X-ray sources. 

For the XMM-Newton observations, we used only the data 
from the pn camera.° We used the Science Analysis System 
(SAS; v1.2) to process the data, following the standard proce- 
dure described in the SAS Data Analysis Threads.’ We used 
the EPPROC tool to produce calibrated event files. A thresh- 
old of 0.4 cts s! was adopted to filter background flares. We 
created good-time-interval files using the TABGTIGEN script, 
and we generated cleaned event files using the EVSELECT 
tool. For each observation, we used the EVSELECT tool to 
extract source and background spectra, with a 30"-radius cir- 
cular source region centered on the optical position of the 
quasar and a 80"-radius circular, source-free background re- 
gion on the same CCD chip as the source region. 

For the 2015 XMM-Newton observation of PHL 1811, 
we also used the Optical Monitor (OM) photometric data 
for constructing its SED. We generated nine exposures for 
the five filters (UVW2, UVM2, UVWI, U, B) using the 
OMICHAIN script, and the photometric measurements of ev- 
ery exposure were recorded in the SWSRLI files. We ex- 
tracted the magnitude measurements from these files and 


5 PHL 1811 has another Chandra observation in 2012 with an exposure time 
of 2.0 ks (see Footnote 18 of Luo et al. 2015); we do not use this short 
Observation in the present study as it does not have sufficient statistics to 


place meaningful constraints on relevant parameters. 


6 We have checked the MOS data, which have lower photon statistics, espe- 
cially in the 5-10 keV band that is of interest to this study; combining the 
pn spectra with the low signal-to-noise ratio MOS spectra might introduce 


additional systematic uncertainties. 


7 http://www.cosmos.esa.int/web/xmm-newton/sas-threads. 


8 The 2004 OM observation of PHL 1811 used only one filter and its photo- 
metric measurement does not suggest any variability compared to the 2015 
OM results. 


adopted the mean magnitude for each filter. The coincidence 
loss corrections (due to the high count rates) were applied 
during the pipeline processing. 


3. X-RAY AND MULTIWAVELENGTH PROPERTIES 


3.1. Soft X-ray Weakness and Obscuration Signatures from 
Archival X-ray Observations 


3.1.1. PG 1001 


PG 1001 was observed by XMM-Newton and Chandra in 
2003 and 2010 with exposure times of 8.6 ks and 4.9 ks, re- 
spectively, and the results were presented in Schartel et al. 
(2005) and Saez etal. (2012). As described in Section 
2.3, we reduced these data and extracted the correspond- 
ing X-ray spectra. We performed simple power-law spec- 
tral fiting of the 2003 XMM-Newton spectrum in the 1— 
10 keV band and the 2010 Chandra spectrum in the 1-8 keV 
band; a lower energy bound of 1 keV was adopted here be- 
cause there is apparent soft X-ray excess emission in the 
0.3-1 keV XMM-Newton spectrum that is probably related 
to ionized absorption (Schartel et al. 2005). The ene 
power-law photon indices are 0.8 +0.2 and 1.3*11, indica- 
tive of X-ray obscuration. From the best-fit results, we de- 
rived two Qox values for PG 1001, and they are shown in 
the aox versus L A plane i in Figure 1; the L, soo A mea 
surement is from pA et al. (2011). Besides the 2 Soft X-ray 
weakness, there is also X-ray flux variability between these 
two observations. We then computed the Aoox parame- 
ters, which are —0.58+0.15 for the XMM-Newton obser- 
vation and —0.75 + 0.16 for the Chandra observation, cor- 
responding to X-ray weakness factors of 32 and 8814 
at rest-frame 2 keV, respectively. The Aoox uncertainties 
were dominated by the oox rms scatter (~ 0.15; Table 5 of 
Steffen et al. 2006) of the OoX-L, o o relation. We note that 
there does not appear to be any strong long-term UV/optical 
variability of the three quasars in our study (maximum flux 
variability amplitudes ~ 16%-80%; see Section 3.3 below), 
and thus the large X-ray weakness factors derived using non- 
simultaneous X-ray and UV/optical data are not significantly 
affected by any UV/optical variability. 


3.1.2. PG 1254 


PG 1254 has a 36.0 ks Chandra observation in 2000 that 
was studied in Sabra & Hamann (2001). We performed 
simple power-law spectral fitting of this Chandra spectrum, 
and the resulting photon index is 0.6+ 0.3, indicative of 
X-ray obscuration. Its location in the aox versus E un A 
plane is shown in Figure 1; the L ne measurement 1s 
from Shen et al. (2011). The corresponding Aoox value is 
—0.74 + 0.16, indicating an X-ray weakness factor of 87:219 
at rest-frame 2 keV. 


3.1.3. PHL 1811 
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Figure 1. X-ray-to-optical power-law slope parameter (ox) vs. 
2500 À monochromatic luminosity for the three quasars, showing 
their significant soft X-ray weakness derived using archival Chan- 
dra and XMM-Newton observations. For PHL 1811, the aox values 
from the 2001 December 5 and December 17 Chandra observations 
are shown; the aox measurements from the two XMM-Newton ob- 
servations not shown here are very close to the 2001 December 5 
Chandra value. The solid blue line represents the aox—L o rela- 


2500 A 
tion from Steffen et al. (2006). The small gray dots and downward 


arrows represent the aox values and upper limits of the Steffen et al. 
(2006) AGN sample, respectively. 


PHL 1811 has two Chandra observations and two 
XMM-Newton observations, as listed in Table 1; the first 
three were presented in Leighly et al. (2007b). From sim- 
ple power-law spectral fitting of the two Chandra and two 
XMM-Newton spectra, we obtained steep spectral shapes 
with photon indices in the range of ~2.0-2.6, consistent with 
those in Leighly et al. (2007b). The two oox values from the 
2001 December 5 and December 17 Chandra observations 
are shown in Figure 1, using the L 0 A value adopted from 
Leighly et al. (2007b). The Aaox values are —0.78 +0.15 for 
the December 5 observation and —0.58 + 0.15 for the Decem- 
ber 17 observation, corresponding to X-ray weakness factors 
of 10815 and 33 at rest-frame 2 keV, respectively. We do 
not show the oox values from the two XMM-Newton obser- 
vations in Figure 1, as they are very close to the 2001 De- 
cember 5 Chandra value (e.g., Aaox = —0.81 for the 2015 
XMM-Newton observation). 

In summary, previous soft X-ray observations of the three 
quasars have revealed significant soft X-ray weakness, with 
the rest-frame 2 keV fluxes ~ 32-108 times weaker com- 
pared to the expectations from their optical/UV emission. 
The soft X-ray spectra of PG 1001 and PG 1254 have flat 
spectral shapes, indicative of X-ray obscuration, while the 
soft X-ray spectra of PHL 1811 do not show evidence for 
X-ray obscuration. 


3.2. Hard X-ray Weakness and Spectral Shape Constraints 
from NuSTAR 


For each NuSTAR observation, we performed aperture 
photometry using the co-added FPMA + FPMB images in the 
three bands (soft, hard, and full). In each image, we extracted 
source counts (S) and background counts (B) from the same 
source and background regions used in the spectral extraction 
in Section 2.2. The encircled-energy fraction of the source 
region is 63.9% according to the NuSTAR point-spread func- 
tion. We determined the source significance by calculating 
the binomial no-source probability (Pp; e.g., Luo et al. 2013), 
which is defined as 


N! 


X(1. p\N-X 
XIND” (L9) ^. (1) 


N 
PX > S)= M7 
X-$ 
In this expression, N = S4 B and p = 1/(1+BACKSCAL), 
where BACKSCAL is the ratio between the exposure-time 
weighted areas of the background and source regions. A 
smaller Pg value indicates a more significant signal. We con- 
sidered the source detected in a band if the measured Pg value 
is smaller than 0.01 (corresponding to a > 2.6c significance 
level). With this criterion, PG 1001 and PG 1254 were not 
detected in the hard band in their 2013 observations, and the 
three quasars were detected in all the other images; these re- 
sults are consistent with the WAVDETECT results in Section 
2.2. In the hard band, the Pg values for PG 1001 and PG 
1254 in their latest observations are 5.4 x 1077 (5.00) and 
3.8 x 1075 (4.60), respectively, and Pg is 1.1 x 10? (3.30) 
for PHL 1811; these values indicate significant detections in 
the hard band. The NuSTAR hard-band images of the there 
quasars are displayed in Figure 2. 

For the detected sources, we computed their 
aperture-corrected net counts (S — B/BACKSCAL)/0.639. 
The associated errors were derived from the 1c Poisson er- 
rors of the extracted source and background counts (Gehrels 
1986). Compared to PG 1001 and PG 1254, PHL 1811 has 
larger relative count errors in the hard band, consistent with 
its larger hard-band Pg value (lower detection significance). 
For undetected sources, we calculated 90% confidence-level 
upper limits on the source counts following the Bayesian 
approach of Kraft et al. (1991). The net counts and upper 
limits in the three bands are listed in Table 2. 

For each quasar in each observation, we derived an effec- 
tive power-law photon index (Teff) from the band ratio, which 
is the ratio between the hard-band (8—24 keV) and soft-band 
(3-8 keV) counts, based on the following procedure. (1) For 
a given set of I values, we produced a set of mock power-law 
spectra using the XSPEC FAKEIT routine and the spectral re- 
sponse files. (2) For each mock spectrum, we computed the 
corresponding band ratio. (3) We interpolated the I’ versus 
band ratio set to derive the I’. value from the measured band 
ratio. The Teg values are listed in Table 2. The lo errors on 
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Figure 2. NuSTAR hard-band (8-24 keV) images of the three 
quasars, smoothed with a Gaussian kernel with standard deviation 
of 2.0 pixels. Each image is centered on the corresponding X-ray 
source position, and the circle indicates the 35"-radius aperture used 
for source extraction. Due to statistical fluctuations, the smoothed 
images of these faint sources do not have circular morphologies fol- 
lowing the NuSTAR on-axis point spread function. 
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Pet were propagated from the errors of the band ratios de- 
rived using BEHR (Park et al. 2006). If the quasar was not 
detected in the hard band, we computed the lower limit on 
Tete from the upper limit on the band ratio calculated using 
BEHR. In this case, [et = 2.0 was adopted in the follow- 
ing calculations of fluxes, flux densities, and luminosities; 
we consider this more appropriate than using the lower limit 
value, and adopting a value different from 2.0 would not af- 
fect the results significantly. 

To compute fluxes, we obtained conversion factors from 
count rates to fluxes in the three bands using the mock spec- 
trum with a photon index of Te. Flux errors were propa- 
gated from the count errors, and flux upper limits were de- 
rived from the upper limits on the net counts. The luminosity 
in the rest-frame 2-10 keV band (Lx) in each observation 
was derived from the full-band flux adopting a power-law 
spectrum with a photon index of Tee. The X-ray fluxes and 
luminosities are listed in Table 2. For PG 1001 and PG 1254 
that have two NuSTAR observations, the photometric prop- 
erties from the two observations are overall consistent within 
the errors, except for the Ief constraints of PG 1001 which 
suggest possible spectral-shape evolution. We stacked their 
photometric measurements to derive average properties; the 
results are also listed in Table 2. 

We calculated the factor of X-ray weakness at rest-frame 
8 keV (fw), which is defined as the ratio between the expected 
and observed 8 keV flux density (fy = f, expected / fu, observed): 
The observed 8 keV flux density was computed from the full- 
band flux for a power-law spectrum with a photon index of 
Petr, and the expected 8 keV flux density was calculated from 
the Steffen et al. (2006) Oox-L 4 A relation adopting a I' = 
2 power-law spectrum. The fw uncertainties were propagated 
from the aox rms scatter of the oox-L ; relation. We 
caution that fw is different from the factor of X-ray weakness 
quantified by the Aoox parameter (Section 1), which is for 
rest-frame 2 keV. 

Compared to the previous hard-band non-detections, the 
deeper NuSTAR observations of PG 1001 and PG 1254 
improve the source-detection significance and provide bet- 
ter constraints on the hard X-ray weakness factors and 
hard X-ray spectral shapes. The NuSTAR observation of 
PHL 1811 provides first-ever hard X-ray constraints for this 
extreme quasar. PG 1001 and PHL 1811 are significantly 
hard X-ray weak (fw + 26-74), while PG 1254 is only X-ray 
weak by a factor of ~ 2.7 in the hard X-rays despite its 
significant soft X-ray weakness (Figure 1). Although the 
2013 observation of PG 1001 suggests a potentially typical 
hard X-ray spectral shape (Tef > 1.5), the 2020 observa- 
tion and the 2013 4 2020 stacked data suggest flat spectral 
shapes (Terr ~ 0.4—1.0). Both observations of PG 1254 sug- 
gest nominal spectral shapes with De ~ 1.8. The PHL 1811 
observation surprisingly reveals that its hard X-ray photon 
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index (Ief = LATS), though loosely constrained, appears 
marginally smaller than its soft X-ray (0.3-5 keV) photon 
index (2.3 + 0.1) from the 2004 XMM-Newton observation 
(Leighly et al. 2007b). 


3.3. Spectral Energy Distributions and Optical/Infrared 
Variability 


We constructed infrared (IR) to X-ray SEDs for the three 
quasars. We collected IR-UV photometric data from the 
Wide-field Infrared Survey Explorer (WISE; Wright et al. 
2010), Near-Earth Object WISE (NEOWISE; Mainzer et al. 
2014), Two Micron All Sky Survey (2MASS; Skrutskie et al. 
2006), Sloan Digital Sky Survey (SDSS; York et al. 2000), 
and Galaxy Evolution Explorer (GALEX; Martin et al. 2005) 
catalogs. For the two PG quasars, we also included their SED 
data from Neugebauer et al. (1987). We added Spitzer photo- 
metric measurements for PG 1001 (Veilleux et al. 2009). For 
PHL 1811, we included the 2001 HST STIS UV spectrum 
and the 2015 XMM-Newton OM measurements. The opti- 
cal and UV data have been corrected for Galactic extinction 
following the de-reddening approach in Cardelli et al. (1989) 
and O'Donnell (1994). The SEDs are shown in Figure 3. 
For comparison, we also plotted in each panel the mean SED 
of high-luminosity radio-quiet quasars in Krawczyk et al. 
(2013) normalized to the 2500 À luminosity. The X-ray com- 
ponent of the mean quasar SED is a [ = 2 power-law contin- 
uum that follows the aox-L À relation. The IR-to-UV 
SEDs of the three quasars are broadly consistent with those 
of typical quasars; the slight deviations in the IR for PG 1001 
and PHL 1811 are within the scatters (~ 0.2—0.25 dex) of the 
mean quasar SED at these frequencies. 

We added soft and hard X-ray measurements to the SEDs. 
We used the 2 keV luminosities from the power-law spectral 
fitting of the Chandra or XMM-Newton spectra (Section 3.1). 
From the NuSTAR full-band fluxes (Section 3.2; stacked re- 
sults were used for PG 1001 and PG 1254), we derived 8 keV 
and 15 keV luminosities adopting power-law spectra with 
the measured Pef values. Compared to the typical quasar 
SED which follows the oox—-L ; relation, the soft and 
hard X-ray weakness of the three quasars is evident. For 
PHL 1811, we also show the soft X-ray spectral slopes from 
the two XMM-Newton observations and the hard X-ray spec- 
tral slope constrained from the NuSTAR observation. The 
spectral slopes differ beyond the 1c level, suggesting that 
X-ray obscuration might also be present in PHL 1811. The 
SEDs of the three quasars also show clearly that they deviate 
from the observed Lx—Lmr relations for typical quasars (e.g., 
Lutz et al. 2004; Mateos et al. 2015; Stern 2015; Chen et al. 
2017; Martocchia et al. 2017), with Lyr being the mid-IR 
luminosity measured at rest-frame 6 um. The offsets from 
the relations are approximately quantified by the fẹ values 
(Table 2), and PHL 1811 is z 70 times X-ray underlumi- 


nous compared to its mid-IR luminosity. Integrating the IR- 
to-X-ray SEDs, we estimate the bolometric luminosities to 
be 3.4 x 10? erg s^! for PG 1001, 2.1 x 10% erg s^! for 
PG 1254, and 4.2 x 10*6 erg s! for PHL 1811. These values 
are consistent with those provided in Shen et al. (2011) and 
Leighly et al. (2007b). 

The IR-to-UV SED data are not simultaneous and may be 
affected by variability. Mild variability is apparent in the 
SED of PG 1001, where the more recent optical and near-IR 
measurements are ~ 2—60% lower than the Neugebauer et al. 
(1987) data. To investigate the optical variability of these 
three quasars, we further examined their long-term opti- 
cal light curves constructed using the public catalogs of 
the Zwicky Transient Facility Data Release 9 (ZTF DR9; 
Bellm et al. 2019) and the Catalina Real-Time Transient Sur- 
vey (CRTS; Drake et al. 2009). The maximum flux variabil- 
ity amplitudes in the ZTF g band range from 13% to 40%, 
and in the CRTS V band they range from 16% to 80%; 
PG 1001 varied the most among the three quasars. Since 
PG 1001 has the strongest optical variability, we also checked 
its IR light curve from the NEOWISE catalog. Between 2014 
May and 2020 November, its maximum flux variability am- 
plitude in the W1 (W2) band is 18% (13%). The mild opti- 
cal/IR variability observed in these quasars suggests that the 
overall accretion power did not change significantly (e.g., by 
factors of > 2) over the years. Compared to the soft and hard 
X-ray weakness factors (A^oox and fw in Sections 3.1 and 
3.2), the UV/optical variability factors are much smaller, in- 
dicating that the X-ray weakness factors assessed using non- 
simultaneous UV/optical data are not heavily biased. 


4. DISCUSSION 


The three quasars in this study were considered among the 
best candidates for intrinsically X-ray weak quasars, with 
coronae that produce weaker X-ray emission than expected 
from the oox-L s relation. Intrinsic X-ray weakness and 
X-ray obscuration are not mutually exclusive. In fact, at least 
for PG 1001 and PG 1254, the presence of X-ray obscuration 
is clearly indicated by their flat spectral shapes in the soft 
X-rays (Section 3.1). Therefore, it was proposed in Luo et al. 
(2014) that the Chandra, XMM-Newton, and NuSTAR data 
of PG 1001 and PG 1254 could be explained by intrinsically 
weak X-ray continua modified by Compton-thin absorption. 
In the subsections below, we first suggest that X-ray obscu- 
ration is also present in PHL 1811 given the hard X-ray con- 
straints. We then propose that, as an alternative to the intrin- 
sic X-ray weakness + X-ray obscuration scenario, the soft 
and hard X-ray weakness of these quasars can be uniformly 
explained under an X-ray obscuration-only scenario, without 
invoking the extra mechanism of intrinsic X-ray weakness. 
The X-ray absorber in this case might be a clumpy accretion- 


Table 2. NuSTAR Photometric Properties 


Object Name Obs. Year Net Source Counts^ Tu 
3-4 3-8 8-24 
keV keV keV 
PG 1001 2013 sate 49h «26 > 1.5 (2.0) 
PG 1001 2020 172334 5753 1133; — 0438 
PG 1001 201342020 22348 9899 125:3 1.0503 
PG 1254 2013 6391 355 «50 > 0.4 (2.0) 
PG 1254 2019 27533 17/38 10432 1.803 
PG 1254 201342019 3384 20653} 13259 1.80.3 


PHL 1811 2015 11323. 5999 ss 1.405 


Flux (10 erg cm? s!) log Lx (erg s!) fs 


3-24 3-8 8-24 2-10 


keV keV keV keV 
7.823 4.413 «5.9 42.7 16525 
7.816 L1 6.6718 42.0 34559 
724+1.2 T3 5.5533 42.2 265 
6.220 2.53 «73 44.5 3.452 
10.04 4345027 57i 44.5 2.5178 
9.12 3.8406 5.3412 44.5 PARE 
7.0511 2:708 4.9115 42.6 74:9 


“The errors were derived from the lo errors of the extracted source and background counts (Gehrels 1986). For undetected sources, we 
calculated 90% confidence-level upper limits on the source counts following the Bayesian approach of Kraft et al. (1991). 


P Effective power-law photon index. If the source is not detected in the 8—24 keV band, a lower limit value is provided, but Ter = 2.0 (as 
shown in parentheses) was adopted in calculating the fluxes, flux densities, and luminosity. 


€ Factor of X-ray weakness at rest-frame 8 keV, derived by comparing the observed 8 keV flux density to that expected from the Steffen et al. 


(2006) aox-L,. s relation assuming a [ = 2 power-law spectrum. The uncertainty is dominated by the scatter of the aox—L 


. 00 A 
relation. 


disk wind that provides variable partial-covering absorption 
to the nuclear X-ray emission. 


4.1. Presence of X-ray Obscuration in PHL 1811 


PHL 1811 was considered an intrinsically X-ray weak 
quasar without any absorption based mainly on the follow- 
ing properties: (1) significantly weak emission with a steep 
spectral shape (I = 2.3 + 0.1) in the 2004 XMM-Newton 
0.3-5 keV spectrum, and (2) flux variability by a factor of 
x~ 4 in 12 days (Leighly et al. 2007b). Hard X-ray (> 5 keV) 
data were not available previously. The 2004 XMM-Newton 
spectrum of PHL 1811 is dominated by background above 
5 keV, with a Pg value (see Section 3.2) of only 0.1 in the 
5-10 keV band. In the deeper 2015 XMM-Newton observa- 
tion, PHL 1811 was significantly detected in the 5-10 keV 
band, with a Pg value of 1.0 x 10^ (3.90), allowing investi- 
gation of its hard X-ray properties. We first derived its pho- 
tometric properties following the same procedure described 
in Section 3.2, though applied to the 0.3-10 keV band in- 
stead of the 3—24 keV band as for the NuSTAR data. The net 
source counts (without aperture correction) in the 0.3-2 keV 
and 2-10 keV bands are 872-35 and IMS respectively, and 
the resulting 0.3-10 keV effective photon index (Tefr0.3-10) 
is 2.0 3: 0.1 from the band ratio of the above two bands and 
the spectral response files. If adopting a different set of en- 
ergy bands, the 0.3-5 keV and 5-10 keV bands, the derived 
counts are 1000+33 and 47*4, and the Tefro.3-10 value be- 


2500 À 


comes 1.8954. The slightly different D'et¢o.3-19 values suggest 


that the spectral shape deviates from a simple power law. 
We then fit the 2015 XMM-Newton 0.3-10 keV spectrum 
with a simple power-law model using XSPEC, under the as- 
sumption that PHL 1811 is an unabsorbed, intrinsically X-ray 
weak source. The spectrum was grouped with at least one 
count per bin, and the W statistic (Footnote 3) was used. 
The resulting photon index is = 2.57 + 0.09; if limiting the 
energy range to 0.3-5 keV, the resulting [ (2.63 + 0.09) is 
consistent within the errors. The data and the best-fit model 
are shown in Figure 4a, and there are significant residuals 
above zz 3 keV. It is possible that a Compton-reflection com- 
ponent from a distant reprocessor (e.g., the torus) contributes 
to the hard X-ray excess emission. We thus added an XSPEC 
pexrav component to the model with its photon index and 
normalization tied to those of the power-law component; the 
only free parameter is the reflection factor (reflection frac- 
tion) and the other parameters are fixed at their default val- 
ues. The best-fit results are shown in Figure 4b. The hard X- 
ray excess can be explained by the Compton-reflection com- 
ponent, but an unrealistically large reflection factor (~ 23) 
is required; the reflection factors for AGN samples are typ- 
ically S 2 (e.g., de Rosa et al. 2008; Ricci et al. 2011, 2017; 
Panagiotou & Walter 2019). We fixed the power-law photon 
index to 2.57 in the above test; when was allowed to vary, 
a larger I (~ 3.0) and an even larger reflection factor (~ 81) 
was derived. Therefore, the 2015 XMM-Newton spectrum of 
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Figure 3. SEDs of the three quasars (see Section 3.3 for the IR- 
UV data). The 2 keV luminosities were derived from the power- 
law spectral fitting of the Chandra or XMM-Newton spectra (Sec- 
tion 3.1). The 8 keV and 15 keV luminosities were derived from 
the NuSTAR photometric results (Section 3.2). For PHL 1811, we 
show the soft X-ray spectral slopes (with negligible uncertainties) 
from the two XMM-Newton observations and the hard X-ray spec- 
tral slope (with the shaded area representing the uncertainties) con- 
strained from the NuSTAR observation, which differ beyond the lo 
level. The green dashed curve in each panel shows the mean quasar 
SED from Krawczyk et al. (2013) normalized to the 2500 A lumi- 
nosity; the X-ray component is a I' = 2 power-law continuum that 


follows the aox—L o relation. 
OX"5500 A 


PHL 1811 cannot be well described by a simple unabsorbed 
power-law continuum plus a reasonable amount of Compton 
reflection. 

In Figure 4a, the hard X-ray residuals peak around 
observed-frame 6 keV (rest-frame z 7 keV), and we thus also 
considered a model where the hard X-ray excess has some 
contribution from a broad Fe Ka line produced via relativis- 
tic disk reflection (e.g., Ross & Fabian 2005; Fabian et al. 
2013). We fit the 2015 XMM-Newton 0.3-10 keV spec- 
trum with the XSPEC re1xill model (Dauser et al. 2014; 
García et al. 2014). The free parameters are I', SMBH spin, 
ionization parameter, inclination angle, and reflection factor 
(reflection fraction), and the other parameters were fixed at 
their default values. This model describes well the spectrum 
with similar residuals to those in Figure 4b, with a best-fit I 
value of 2.02 + 0.06 and a best-fit reflection factor of 10465, 
The I value is smaller than 2.57 because the soft X-rays are 
now dominated by ionized-disk reflection instead of an in- 
trinsic power-law continuum. The reflection factor is still 
large, but unlike the Figure 4b modeling with a distant reflec- 
tor, relativistic reflection from the inner accretion disk could 
produce an extremely large reflection factor if much of the 
coronal emission cannot reach the observer due to the light 
bending effects near the SMBH (e.g., Dauser et al. 2014). 
Nevertheless, a reflection-dominated spectrum with strong 
light bending effects is still in contrast with the scenario of in- 
trinsic X-ray weakness where the observed spectrum should 
be dominated by the intrinsic power-law continuum from a 
weak corona. 

The 3-24 keV spectral shape constrained from the 2015 
simultaneous NuSTAR observation, Deg = [405 (Sec- 
tion 3.2), provides additional support that the spectral shape 
likely deviates from a steep unabsorbed power law. We thus 
applied the above zpow + pexrav model to jointly fit the 
2015 XMM-Newton 4 NuSTAR spectra and investigate if 
the spectra can be explained by an unabsorbed, intrinsically 
weak power-law continuum plus typical Compton reflection. 
The results are shown in Figure 5a, and they are consistent 
with those for the XMM-Newton spectrum alone (Figure 4b), 
requiring a huge reflection factor. We then tested replacing 
the pexrav component with the self-consistent Compton- 
reflection model borus02 (Balokovié et al. 2018) with the 
photon index and normalization tied to those of the power- 
law component. The complete XSPEC model is phabs x 
(zpow+atable{borus02.fits}), and the free param- 
eters are the power-law normalization, photon index (T`), and 
column density (Ny) of the reprocessed component.’ The 


? The other parameters of the borus02 model are fixed at the default values 
including a high-energy cutoff of 300 keV, an inclination angle of 60°, a 
torus covering factor of 0.5 (corresponding to a half-opening angle of 60°), 
and an iron relative abundance of 1. 
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Figure 4. The 2015 XMM-Newton spectrum for PHL 1811 overlaid with the best-fit (a) simple power-law model and (b) simple power-law 
plus Compton-reflection model (with a fixed photon index). The spectra are grouped for display purposes only. The bottom panels display 
the ratios of the data to the best-fit models. Compared to the simple power-law model, there are significant fitting residuals above & 3 keV, 
which require an unrealistically large reflection factor ( 23) to explain. The spectrum thus cannot be well described by a simple unabsorbed 
power-law continuum plus a reasonable amount of Compton reflection. 
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Figure 5. Similar to Figure 4, displaying the 2015 XMM-Newton and NuSTAR spectra for PHL 1811 overlaid with the best-fit power-law 
plus Compton-reflection model. The Compton-reflection component was modeled with (a) pexrav or (b) borus02. An unrealistically large 
reflection factor (~ 24) is required in the pexrav model, while the self-consistent borus02 Compton reflection is not sufficient to account 


for the excess emission in the hard X-rays. 


best-fit results are displayed in Figure 5b. There are still sig- 
nificant fitting residuals above z 5 keV, demonstrating again 
that a typical level of Compton reflection cannot account for 
the excess emission in the hard X-rays. 

The 2015 XMM-Newton and NuSTAR observations of 
PHL 1811 thus reveal that, in addition to a steep power-law 
(T = 2.63 + 0.09) continuum in the 0.3—5 keV band, there is 
significant hard X-ray (= 5 keV) excess emission. The ex- 
cess emission can be modeled with Compton reflection of 
a soft X-ray continuum that is much stronger (by a factor 
of = 20) than the observed one. We note that, from the 
best-fit zpow + pexrav results in Figure 5a, PHL 1811 


was intrinsically X-ray weak by a factor of ~ 140 compared 
to the Steffen et al. (2006) aox-L ; relation. If we in- 
stead assume that PHL 1811 was intrinsically X-ray nor- 
mal (i.e., raising the power-law normalization by a factor 
of 140) and add a scaling factor to the power-law compo- 
nent c * zpow + pexrav, the spectra can be equally well 
described with c = 0.796 (1/140) and a reasonable reflec- 
tion factor of 0.17 (24/140). A physical interpretation of the 
above model is that PHL 1811 was affected by Compton- 
thick absorption and the observed XMM-Newton and NuS- 
TAR spectra were dominated by a fraction of the intrinsic 
continuum scattered by a large-scale highly ionized “mir- 


12 


ror” (fecatter; typically within a few percent; e.g., Turner et al. 
1997; Cappi et al. 2006; Ueda et al. 2007; Winter et al. 2009; 
Yamada et al. 2020; Gupta et al. 2021) in the soft X-rays 
and a reprocessed component from the absorber in the hard 
X-rays. The steep soft continuum could also represent a 
fraction (fieak) of the intrinsic continuum leaking through a 
clumpy absorber. Therefore, instead of employing intrinsic 
X-ray weakness plus an unrealistically large reflection fac- 
tor, we could interpret the 2015 XMM-Newton and NuSTAR 
data of PHL 1811 with Compton-thick obscuration. 

In addition, the second property mentioned above, the 
short-term soft X-ray variability between the two Chandra 
observations in 2001 (e.g., see Figure 1), was used to ar- 
gue against a large-scale scattered component (fécatter above) 
dominating the 0.3-5 keV X-ray spectrum. However, from 
recent investigations of extreme X-ray weakness and ex- 
treme X-ray variability among super-Eddington accreting 
AGNs (e.g., SDSS JO75101.42+291419.1 and samples from 
Liu et al. 2019, 2021; SDSS J081456.10-- 532533.5 from 
Huang J. et al. in prep.), it appears plausible that such 
weak, steep, and variable soft X-ray emission might originate 
from variable fractions (fieak) of leaked intrinsic continuum 
through a large solid-angle, high column-density, clumpy 
(partial-covering) absorber. Thus fast variability and steep 
spectral shapes in the soft X-rays do not necessarily rule out 
X-ray obscuration. In summary, unlike the model considered 
previously where PHL 1811 lacks absorption, hard X-ray 
data suggest that X-ray obscuration may well be present. 


4.2. An Obscuration-Only Scenario Without Intrinsic X-ray 
Weakness 


The 2013 NuSTAR observations of PG 1001 and PG 1254, 
with no hard-band detections or Ief measurements, were 
not sufficiently constraining to establish that intrinsic X-ray 
weakness must be present, which motivated the present study 
with deeper NuSTAR observations. The deeper NuSTAR ob- 
servations now provide hard-band detections and Te mea- 
surements (albeit with large uncertainties). The PG 1001 
spectral shape appeared flatter in the 2020 observation (Terr = 
0.4436 compared to Terr > 1.5), suggesting the presence of 
absorption and probably spectral-shape evolution between 
these two epochs. The stacked [De value (1.0+8:5) from the 
two NuSTAR observations is still small compared to the typ- 
ical value of ~ 2 for an unabsorbed spectrum. The PG 1254 
spectral slope appears typical for unobscured quasars (Teff = 
1.8223 in the 2019 observation), but it is also clear that this 
quasar is X-ray weak by only a factor of a few in the hard 
X-rays (see Figure 3). Considering its significant soft X-ray 
weakness and the flat spectral shape in the soft X-rays (Sec- 
tion 3.1), the nominal hard X-ray spectral shape in the NuS- 
TAR data might be explained by Compton-thin absorption. 
Since PG 1254 is at z = 1.026, we are likely observing the 


penetrating hard X-rays with NuSTAR through an absorber 
with a large but Compton-thin Ny value; this would explain 
the strong X-ray weakness at 2 keV and much reduced X-ray 
weakness at higher energies. As discussed in Section 4.1 
above, hard X-ray data suggest that X-ray obscuration is also 
present in PHL 1811. 

Motivated by these NuSTAR results, we argue that intrin- 
sic X-ray weakness is probably not required to explain the 
extreme X-ray weakness of these quasars. We explore the 
possibility of interpreting universally the NuSTAR, Chandra, 
and XMM-Newton spectra with an obscuration-only scenario 
where these quasars are intrinsically X-ray normal (following 
the aox—L s relation). We investigate below, via XSPEC 
spectral fitting, 1f the multi-epoch soft and hard X-ray spectra 
can be described by nominal-strength X-ray emission modi- 
fied by our adopted obscuration model with the absorber pa- 
rameters (the column density and partial-covering fraction) 
allowed to vary. Since the X-ray data do not have sufficient 
statistics for complex spectral fitting, we had to simplify the 
model and fix many of the model parameters. Moreover, we 
actually cannot rule out the scenario of obscuration - intrin- 
sic X-ray weakness, which has an extra degree of freedom 
(i.e., normalization of the X-ray continuum) compared to the 
obscuration-only scenario. Our focus here is to investigate 
if we can explain the observed X-ray emission without in- 
volving intrinsic X-ray weakness from, e.g., an anomalous 
corona. 

Although the pexrav model appears able to describe the 
PHL 1811 spectra well (Section 4.1 and Figure 5a), it does 
not provide constraints on the absorption column densities. 
We thus employed the self-consistent borus02 model to 
describe the reprocessed component from the absorber. The 
XSPEC spectral model is 


phabs * (zphabs * cabs * C] * zpow 


+ co * zpowt+ atable{borus02.fits}). 


In this model, phabs accounts for the Galactic absorption, 
and zpow is the intrinsic power-law continuum that is X-ray 
normal with respect to the Steffen et al. (2006) aox-L o A 
relation. A large fraction (c1) of the intrinsic continuum is 
modified by heavy neutral absorption (zphabs) and Comp- 
ton scattering (cabs). The absorber is clumpy, allowing a 
fraction (fieak 2c? = 1—c1) of the intrinsic continuum to leak 
through. There is probably also a large-scale scattered com- 
ponent (fiue). Since the leaked component usually dom- 
inates, we treated them together and do not separate fscatter 
from fik in the following study. The reprocessed com- 
ponent from the absorber is modeled with borus02 with 
the normalization and photon index tied to those of zpow. 
We fixed the inclination angle to 60? and allowed the ab- 
sorber covering factor to vary. The other borus02 param- 
eters were fixed at the default values (Footnote 9). We also 


tied the absorption column densities (Ny) in the zphabs, 
cabs, and borus02 components. The emergent spectrum 
is thus a combination of the transmitted (absorbed) com- 
ponent, leaked component, and reprocessed component. A 
schematic illustration of the setup is shown in Figure 6a, and 
an example of the X-ray spectral components for a given set 
of the absorber parameters is displayed in Figure 6b. 

There are a couple of caveats regarding the above model. 
First, the absorber is probably partially ionized instead of 
being neutral. Ionized absorption produces distinctive spec- 
tral features below z 1 keV, but the continuum shapes above 
~ 1 keV are similar to those from neutral absorption unless 
the ionization level is high (e.g., Figure 1 of Netzer 1993). 
For the three quasars studied here, PG 1254 has few photons 
below 1 keV, PHL 1811 shows no absorption signatures be- 
low + 5 keV, and PG 1001 has clear soft X-ray excess emis- 
sion in the 0.3-1 keV band which was interpreted with ion- 
ized absorption (Schartel et al. 2005). Their > 1 keV spec- 
tra do not have sufficient photon statistics to distinguish ion- 
ized absorption from neutral absorption or constrain reliably 
ionization parameters. Therefore, in the above model, we 
adopted neutral absorption for simplicity, and we did not use 
the < 1 keV XMM-Newton or Chandra data for PG 1001. 
The soft excess of PG 1001 in the obscuration scenario will 
be discussed in Section 4.2.1 below. Second, since there is no 
optical/UV extinction, the absorber should not be the torus 
described in the borus02 model with a toroidal geome- 
try. Instead, it is likely a small-scale clumpy dust-free wind 
launched from the accretion disk (see Section 4.3 below). 
Therefore, the borus02 model does not provide an accurate 
description of the reprocessed emission (both the continuum 
and the Fe Ka line) from the absorber (e.g., wind). However, 
since our purpose here is not to recover precise absorber pa- 
rameters but to simply investigate if the obscuration-only sce- 
nario is a valid alternative to the scenario of intrinsic X-ray 
weakness + obscuration, and the current simplified model ap- 
pears able to explain reasonably well the multi-epoch X-ray 
spectra of the three quasars (as discussed below), we defer 
detailed modeling to future studies which likely will require 
much better spectral quality. 

We applied the above model to explain the multi-epoch 
X-ray spectra of the three quasars. We jointly fit the NuS- 
TAR, Chandra, and XMM-Newton spectra for each of the 
three quasars. The free parameters are I', fleak (C2), Nu, 
and the absorber covering factor (cos0,, in Figure 6a); fieak 
and Ny were allowed to vary between the observations while 
the other two parameters were tied. In addition, for the 
latest NuSTAR observation of PG 1001 and the two NuS- 
TAR observations of PG 1254, we tied the ficak parameter 
to that of the non-simultaneous Chandra observation, as the 
NuSTAR spectra are not sensitive to this parameter (high- 
energy spectra do not have a significant leaked component). 
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These quasars are considered to be intrinsically X-ray nor- 
mal, with the intrinsic foxey values fixed at those expected 
from the Steffen et al. (2006) aox-L 9 À relation. The 
fxev and I values define the intrinsic continua (normaliza- 
tions of zpow). For each object, the zpow normalization 
was first fixed at the value derived with I = 2, and then an 
iterative procedure was performed to adjust the zpow nor- 
malization according to the best-fit [ value. A few iterations 
were needed until the intrinsic continuum converged. 

The best-fit parameters are listed in Table 3 and the best-fit 
models are displayed in Figure 7. We also list in Table 3 the 
main components (transmitted, reprocessed, and/or leaked) 
that dominate the emergent spectrum of each observation. 
The I values for PG 1001 and PHL 1811 pegged at 2.6, the 
upper bound allowed by the borus02 model. Several Ny 
values do not have upper errors as the emergent spectrum 
is dominated by the leaked component which is not sensi- 
tive to Nu; a few of these pegged at log Nu /cm ? = 25.5, the 
upper bound allowed by the borus02 model. The fitting 
results are acceptable overall considering the fitting statistics 
(W /dof; Figure 7) and residuals, indicating that the obscu- 
ration scenario is able to explain the multi-epoch X-ray data 
without involving intrinsic X-ray weakness. From the best-fit 
models, we computed the fy values for the NuSTAR spectra, 
and these are also listed in Table 3. These hard X-ray weak- 
ness factors are comparable to those listed in Table 2, and 
any discrepancy is mainly due to the different models used 
(a simple power-law model is assumed in deriving the photo- 
metric results in Table 2). 

Given the limited X-ray data quality and the simplifica- 
tions/assumptions made during the modeling process, we do 
not consider these spectral fitting results fully accurate de- 
scriptions of the absorber properties. Nevertheless, they pro- 
vide important clues for explaining the unusual X-ray prop- 
erties under our proposed obscuration-only scenario. The 
best-fit results are consistent with our qualitative expectation 
above. The multi-epoch spectra of PG 1001 are explained by 
heavy or even Compton-thick obscuration. A strong leaked 
component is required to explain its 2013 NuSTAR spectrum 
as the spectral shape is likely steeper than an absorbed power 
law (Legg > 1.5), while the 2020 NuSTAR spectrum is af- 
fected by typical Compton-thick obscuration. The Chandra 
spectrum of PG 1254 requires very Compton-thick obscu- 
ration due to the significant weakness of this z = 1 spec- 
trum. The NuSTAR observations of PG 1254 are explained 
by heavy but Compton-thin obscuration. PHL 1811 was al- 
most always affected by Compton-thick obscuration, and the 
emergent spectra are largely dominated by the leaked com- 
ponent (though with small fieax values). The reprocessed 
component does not contribute much as the absorber appears 
to have a large covering factor (cos, œ% 0.6) which blocks 
direct reflected radiation from the opposite side of the ab- 
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Figure 6. (a): A schematic diagram of the obscuration scenario. The black dot represents the SMBH, surrounded by the X-ray corona shown 
in pink. The yellow region represents the accretion disk. The green arrows represent the absorber (e.g., an outflowing disk wind) with a half- 
opening angle of Ooa and a covering factor of cos0oa. The system inclination angle is inc. The dashed line traces the transmitted (absorbed) 
radiation through the absorber in the direction of the observer (transmitted component). The dotted curves trace the radiation reflected from 
the opposite side of the absorber and the radiation scattered to the line of sight through the absorber (reprocessed component). The absorber is 
clumpy allowing a fraction (fica) of the intrinsic continuum to leak through without being absorbed (leaked component). There is probably also 
a fraction (fiue) of the intrinsic continuum scattered to the line of sight by a large-scale highly ionized “mirror” (not shown in the cartoon). 
We treated ficak and fecatter together in this model. (b): An example of the X-ray spectral components for a given set of the absorber parameters. 
The black solid curve shows the total observed model spectrum, while the gray long-dashed line, yellow dashed-dotted curve, red dotted curve, 
and blue dashed line represent the intrinsic power law, transmitted component, reprocessed component, and leaked component, respectively. 


sorber (e.g., see the top dotted curve in Figure 6). Therefore, 
PHL 1811 can reach a very large hard X-ray weakness factor 
(fw) in the NuSTAR observation. For the 2015 simultane- 
ous XMM-Newton and NuSTAR observations of PHL 1811, 
the derived parameters differ slightly. We also tried to tie all 
the parameters in these two observations, and the results are 
listed in the last row of Table 3. We do not consider the small 
discrepancy a serious issue as there may be cross-calibration 
uncertainties between the XMM-Newton and NuSTAR data 
(e.g., Madsen et al. 2021). 


4.2.1. Soft X-ray Excess of PG 1001 in the Obscuration-Only 
Scenario 


Soft X-ray excess emission (typically below œ% 1 keV) is 
observed in a large fraction of type 1 AGNs, the origin of 
which is still under debate and may be attributed to ion- 
ized absorption (e.g., George et al. 1998; Gierlifski & Done 
2004), ionized disk reflection (e.g., Ross & Fabian 2005; 
Crummy et al. 2006), or Comptonization in a warm corona 
(e.g., Done et al. 2012). Of the three quasars in this study, 
only PG 1001 shows a clear soft-excess component in both 
its XMM-Newton and Chandra spectra. We thus did not con- 
sider the soft excess in the above modeling, excluding the 
the « 1 keV data for PG 1001. We explore here if the soft 
X-ray excess emission of PG 1001 can be explained in the 
obscuration-only scenario. We focus on the XMM-Newton 
spectrum in the following discussion, as the Chandra spec- 
trum has only 19 counts in the 0.3-8 keV band. Neverthe- 
less, we verified that the Chandra spectrum yields consistent 
results. 


The soft-excess emission of PG 1001 is at a comparable 
flux level to the > 1 keV power-law component that is sig- 
nificantly weak compared to the expectation from the aox— 
Soin A relation (Section 3.1.1). Therefore, the soft excess 
is also significantly weak compared to typical levels. In the 
obscuration-only scenario, PG 1001 has a nominal-strength 
hot corona, and likely also a nominal-strength warm corona. 
A natural interpretation would then be that the soft-excess 
emission is also filtered by the absorber if it is from the warm 
corona. We thus fitted the 0.3-10 keV XMM-Newton spec- 
trum with the same model described in Section 4.2 above 
plus an additional component to describe the soft excess from 
the warm corona. We tested simple power-law (zpow), disk 
multi-black body (diskbb), or Comptonization (compTT) 
models for this soft-excess component, and the three choices 
were all able to describe the spectrum well with comparable 
statistics. In Figure 8, we show the best-fit results with the 
power-law model. The soft excess has a large photon index 
(~ 5.0) and a small normalization (~ 0.7% of the normaliza- 
tion for the intrinsic > 1 keV power-law continuum), and the 
other free parameters (fieak and Ny) are consistent with those 
in Table 3 (first row) within the errors. One interpretation is 
thus that the observed soft excess is the leaked portion of the 
warm corona emission through the same dust-free absorber, 
and the leaked fraction is similar to or even the same as that 
for the main component. 

The soft X-ray excess emission of PG 1001 has also been 
suggested to be due to ionized absorption (Schartel et al. 
2005). We verified that the 0.3-10 keV XMM-Newton 
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Table 3. Best-fit Parameters for the Multi-epoch X-ray Spectral Fitting 


Object Name Observatory Obs. Date COS Üoa 
PG 1001 XMM-Newton 2003 May4 0.54006 
PG 1001 Chandra 2010 Jan 11 - 
PG 1001 NuSTAR 2013 Jun 28 i 
PG 1001 NuSTAR 2020 May 23 - 
PG 1254 Chandra 2000 May 29 0.38016 
PG 1254 NuSTAR 2013 Jun 8 = 
PG 1254 NuSTAR 2019 Jun 8 - 
PHL 1811 Chandra 2001Dec5 0.5800) 
PHL 1811 Chandra 2001 Dec 17 - 
PHL 1811 XMM-Newton 2004 Nov 1 - 
PHL 1811 XMM-Newton 2015 Nov 29 - 
PHL 1811 NuSTAR 2015 Nov 28 E 
PHL 1811 XMM-Newton + NuSTAR^ 2015 Nov - 


Partial-Covering Absorption fi Main Com.? 
fieak T log Ny (cm?) 

1.903 x 10? 2.6o9^ 23.5+0.1 = leak, tra 
7.538 x 10° - 23.75 - leak, tra 
0.10 + 0.05 - 24.305 16.7 leak, rep 
7.5 x 10? (tied) - 24.240.1 52.7 tra, rep 
TART «104 . 242705 25.5-10 - leak, rep 
— (tied) - 24.0+ 0.2 8.7 tra, rep 

- (tied) - 23.602 3.2 tra 
6.005 x 10° 2.6-0.05 24.103 - leak, rep 

(2.82: 0.2) x 10? - 24.9.95 - leak 

(6.7 2: 0.4) x 10? - 25.5.95 - leak 
(7.4 3: 0.3) x 10? - 24.6.02 z leak, rep 
(1.1 +0.8) x 10? - 24.8+0.2 98.3 leak, rep 
(7.4 0.3) x 10? - 24.804 - leak, rep 


“The factor of X-ray weakness at rest-frame 8 keV derived from the best-fit model, for comparison with the results in Table 2. 


P The “Main Com.” column lists the dominant component/components in the emergent spectrum: “leak” represents the leaked/scattered compo- 
nent, "tra" represents the transmitted (absorbed) component through the absorber (dashed line in Figure 6), and “rep” represents the reprocessed 


component from the absorber (dotted curves in Figure 6). 


Tn this case, the fitting parameters for the XMM-Newton and NuSTAR observations were tied. 


d A value without an upper error is effectively bound by the allowed upper bound of the borus02 model. 


spetrum can be acceptably fitted with a simple partial- 
covering ionized absorption model (zxipcf*zpow), fix- 
ing I’ = 2.6 and the power-law normalization at the X-ray 


nominal value from the Oox-L, 9 A relation. The resulting 


ionization parameter is € ~ 91 erg cm s! with Ny c 4.7 x 
10? cm"? and a covering fraction of ~ 99.1%. Replacing the 
neutral absorption (zphabs) in the Section 4.2 model with 
zxipcf yields consistent results, as the reprocessed com- 
ponent (borus02) is not important in the XMM-Newton 
spectrum. The soft excess can thus also be explained with 
ionized absorption, which is possible considering that the ab- 
sorber (e.g., disk wind) is probably partially ionized. Overall, 
we consider that our proposed obscuration-only scenario can 
plausibly explain the soft X-ray excess emission of PG 1001. 


4.3. Clumpy Accretion-Disk Wind as the Absorber and 
Implications 


Under our proposed obscuration scenario above, PG 1001 
and PG 1254 are BAL quasars with intrinsically normal 
X-ray emission. Thus, they are probably similar to the typi- 
cal BAL quasars that generally show X-ray obscuration. Be- 
sides the few low-redshift BAL quasars with NuSTAR obser- 


vations (Luo et al. 2013, 2014; Teng et al. 2014), a few high- 
redshift BAL quasars have been suggested to be intrinsically 
X-ray weak (Liu et al. 2018). From the systematic Chan- 
dra survey of the 29 high-ionization LBQS BAL quasars at 
z ^: 2, two intrinsically X-ray weak candidates were identi- 
fied based on their nominal spectral shapes (Teg ~ 2) and 
significant hard X-ray weakness factors (fy ~ 12-15); at 
these redshifts, Chandra observations were able to provide 
rest-frame hard X-ray constraints over a similar band to the 
NuSTAR observations of low-redshift objects. We recently 
obtained a long XMM-Newton observation for one of the 
candidates, LBQS 1442 — 0011. The observation was heav- 
ily affected by background flares and the cleaned exposure 
time is only 4096 of the total. The results are summarized 
in the Appendix. The XMM-Newton observation suggests 
that the spectral shape became flatter (Tef ~ 1) and the fac- 
tor of hard X-ray weakness (fy) also dropped to ~ 4. We 
thus consider that the X-ray weakness of this quasar can 
also be described by our proposed scenario of variable ob- 
scuration, without invoking intrinsic X-ray weakness. It is 
natural to consider that the X-ray absorbers in PG 1001, 
PG 1254, and LBQS 1442 —0011 are similar to those in typi- 
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Figure 7. The Chandra, XMM-Newton, and NuSTAR spectra for the three quasars, overlaid with the best-fit models using the obscuration 
model described in Section 4.2. The spectra are grouped for display purposes only. The bottom panels display the ratios of the data to the 
best-fit models. We break the five observations of PHL 1811 into two panels for clarity. For PG 1001, the energy ranges used are 1-8 keV for 
the Chandra observation and 1-10 keV for the XMM-Newton observation due to the apparent soft X-ray excess emission at lower energies. 
The obscuration scenario explains reasonably well the multi-epoch broad-band X-ray spectra of the three quasars. 


cal BAL quasars; i.e., the shielding gas or clumpy accretion- 
disk wind. They might have an extreme version of the ab- 
sorber in terms of its high column density and/or large cov- 
ering factor. 

An obscuration explanation for the X-ray weakness of 
PHL 1811 would connect it with the PHL 1811 analogs 
studied in Luoetal. (2015). PHL 1811 and its analogs 
belong to a broader category of quasars, weak emission- 
line quasars (WLQs), a small population of type 1 quasars 
that show unusually weak UV emission lines (e.g., C Iv). 
Systematic X-ray surveys of WLQ samples have revealed 
that a large fraction (~ 3096-5096) of them are X-ray 
weak (e.g., Luo et al. 2015; Ni et al. 2018, 2022; Pu et al. 
2020). These WLQs have typically been selected to lie 
at z = 1.5, and thus Chandra or XMM-Newton observa- 
tions provide rest-frame hard X-ray constraints. The in- 
dividual and stacked effective power-law photon indices 
for the X-ray weak WLQs are in general flat (Tef ~ 1.2), 


suggesting an obscuration scenario. The absorber is pro- 
posed to be the geometrically thick inner accretion disk 
and/or its associated disk wind, which also shields the nu- 
clear EUV/X-ray ionizing radiation from reaching the BELR, 
causing the weak emission lines. Thick inner accretion 
disks are expected in WLQs as they are considered to have 
high or even super-Eddington accretion rates that result in 
thick disks (e.g., Abramowicz et al. 1988; Mineshige et al. 
2000; Wang & Netzer 2003; Ohsuga & Mineshige 2011; 
Jiang et al. 2014). Powerful accretion-disk winds launched 
via radiation pressure are expected in such systems (e.g., 
Takeuchi et al. 2014; Jiang etal. 2019; Giustini & Proga 
2019). PHL 1811 likely has a super-Eddington accretion rate, 
and the large covering-factor absorber is probably the strong 
wind associated with a thick accretion disk. 

Therefore, the three quasars in this study likely share 
a similar nature, with partial-covering absorption from 
clumpy dust-free winds (e.g., Figure 6). X-ray absorp- 
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Figure 8. The 0.3-10 keV XMM-Newton spectrum of PG 1001 
fitted with the Section 4.2 model plus an additional power-law com- 
ponent. The spectrum is grouped for display purposes only. The 
bottom panel displays the ratios of the data to the best-fit model. 
The additional power-law component describes well the soft-excess 
emission with a large photon index of z 5.0. The soft excess is 
probably the leaked portion of the warm corona emission through 
the same dust-free absorber that obscures the > 1 keV main compo- 
nent. 


tion from clumpy winds/outflows has been observed in typi- 
cal type 1 AGNs (e.g., Kaastra et al. 2014; Mehdipour et al. 
2017, 2021; Dehghanian et al. 2019; Laha et al. 2021), al- 
though the absorption strength is often not comparable to 
the extreme X-ray weakness found in our three quasars. 
The wind strength and density likely have an Eddington ra- 
tio dependence, as heavy or even Compton-thick absorp- 
tion has been observed in super-Eddington accreting AGNs 
(e.g., Longinotti et al. 2019; Liuetal. 2021). The three 
quasars in this study probably have super-Eddington accre- 
tion rates that drive powerful and high-density clumpy winds. 
They have large estimated Eddington ratios as listed in Ta- 
ble 1. We note that these Eddington ratios do not repre- 
sent accurately the accretion power in the super-Eddington 
regime as a large fraction of the power may be advected 
into the SMBH or converted into mechanical energy of the 
wind (e.g., Jiang et al. 2019); also see Section 2.1 for fur- 
ther discussion. The large intrinsic X-ray photon indices de- 
rived from spectral fitting (T ~ 2.1-2.6; Table 3) also sug- 
gest super-Eddington accretion rates (e.g., Shemmer et al. 
2008; Huang et al. 2020). The NLQI classification and the 
weak [O mI] emission of PG 1001 and PHL 1811 (Sec- 
tion 2.1) provide additional support of super-Eddington ac- 
cretion in these two quasars (e.g., Boroson & Green 1992; 
Sulentic et al. 2000; Shen & Ho 2014). 

It is somewhat odd that PG 1001 shows a strong C Iv 
emission line in its UV spectrum (Brandt et al. 2000), as the 
strong and high-density wind and/or the thick inner accretion 
disk should be able to shield the BELR from the nuclear ion- 
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ization, resulting in a WLQ like PHL 1811. It also appears 
unusual that PHL 1811 does not show any significant UV ab- 
sorption lines (i.e., BALs). Perhaps the dynamical nature of 
the wind (e.g., variable Ny and covering factor) causes the 
apparent discrepancy, and multi-epoch UV spectra might be 
able to shed some light. For example, a z ~ 2 WLQ has re- 
cently been found to undergo BAL transformation (Yi et al. 
2022). Geometric effects might also play a role, as the line of 
sight to the X-ray corona, line of sight to the accretion-disk 
UV continuum region, and the direction from the nucleus to 
the BELR are different from each other, and thus the emer- 
gent X-ray and UV spectra depend on the physical configu- 
ration of the clumpy wind (e.g., Giustini & Proga 2019). 

Although our proposed obscuration scenario was based 
on the new sensitive NuSTAR observations of these three 
quasars, the general connection with obscuration from the 
disk wind suggests that this scenario may be applicable to the 
other intrinsically X-ray weak quasar candidates (e.g., those 
in Nardini et al. 2019 and Laurenti et al. 2022). Obscuration 
from the clumpy disk wind would predict X-ray variability 
from varying obscuration. For the three quasars in this study, 
PG 1001 and PHL 1811 showed clear soft X-ray variability 
(e.g., Figure 1); the 12-day variability timescale of PHL 1811 
does not provide any strong constraints on the wind velocity 
as a wind clump only needs to move a fraction of the corona 
size. The NuSTAR observations of PG 1001 suggest some 
hard X-ray variability, at least in the spectral shape (Table 2). 
The PG 1254 NuSTAR observations do not provide sufficient 
photon statistics to identify hard X-ray variability. In addi- 
tion, LBQS 1442—0011 likely has hard X-ray variability (see 
the Appendix). The WLQs have limited multi-epoch obser- 
vations, and a few of them have been found to vary strongly 
between X-ray normal and X-ray weak states (Miniutti et al. 
2012; Nietal. 2020). A small fraction of super-Eddington 
accreting AGNs have also been found to vary between X- 
ray normal and X-ray weak states (e.g., Liu et al. 2019; 
Boller et al. 2021; Liu etal. 2021), and a few also show 
steep spectra in the low state. Another characteristic prop- 
erty of such X-ray variability is that there is no contempo- 
raneous optical/UV continuum or emission-line variability, 
which argues against changes of accretion rates and supports 
the obscuration scenario. This property also makes these 
AGNS distinct from the unusual population of “changing- 
look" AGNSs (e.g., LES 1927 +654; Trakhtenbrot et al. 2019; 
Ricci et al. 2021) that also show extreme X-ray variability 
but are generally attributed to changes of accretion rates or 
tidal-disruption events. Multi-epoch X-ray observations of 
the intrinsically X-ray weak quasar candidates might be able 
to reveal X-ray variability and help clarify their nature. 


5. SUMMARY AND FUTURE WORK 
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In this paper, we used NuSTAR observations of PG 1001, 
PG 1254, and PHL 1811 to constrain their hard X-ray 
(= 5 keV) weakness and spectral shapes, and thus to inves- 
tigate the nature of their extreme X-ray weakness. These 
quasars show very weak soft X-ray emission (Figure 1), and 
they were previously proposed to be intrinsically X-ray weak, 
with the X-ray coronae producing weak continuum emission 
relative to their optical/UV emission (deviating below the 
aox-L, . "1 relation). The multi-epoch soft and hard X-ray 
observations are summarized in Table 1. NuSTAR aperture 
photometry was presented in Section 3.2, and the results are 
summarized in Table 2 and Figure 3. The NuSTAR spectral 
shapes for PG 1001 and PHL 1811 appear flat (Ce = 10:07 
and Teg = L5. respectively), while the shape is nominal 
for PG 1254 (Tere = 1.80.3). PG 1001 and PHL 1811 
are significantly hard X-ray weak compared to the expecta- 
tions from their optical/UV emission (fy at 8 keV ~ 26-74), 
while PG 1254 is only X-ray weak by a factor of ~ 3. The 
PHL 1811 hard X-ray photon index appears smaller than its 
soft X-ray (0.3—5 keV) photon index (2.3 30.1). Spectral 
modeling suggests that its 2015 XMM-Newton and NuSTAR 
spectra cannot be described by an intrinsically weak contin- 
uum plus a reasonable amount of Compton reflection (Sec- 
tion 4.1 and Figures 4 & 5). 

In light of the new NuSTAR results, a variable X-ray ab- 
sorber can account for all the observations of these X-ray 
weak quasars. We propose that, as an alternative to the in- 
trinsic X-ray weakness + X-ray obscuration scenario, the soft 
and hard X-ray weakness of these quasars can be uniformly 
explained under an X-ray obscuration-only scenario, with- 
out invoking the extra mechanism of intrinsic X-ray weak- 
ness (Section 4.2). In this scenario, the weak emergent spec- 
trum is a combination of the transmitted component modi- 
fied by absorption, the leaked component through a clumpy 
absorber (including a distant scattered component), and the 
reprocessed component reflected/scattered from the absorber 
(Figure 6). This partial-covering absorption scenario pro- 
vides adequate explanations of the multi-epoch X-ray data 
of these quasars, and the X-ray variability is mainly induced 
by the varying column density and leaked fraction (partial- 
covering fraction) of the absorber (Table 3 and Figure 7). 

We propose that the absorber is a clumpy dust-free wind 
launched from the accretion disk (Section 4.3). These 
quasars probably have super-Eddington accretion rates which 
result in geometrically thick inner accretion disks and pow- 
erful winds with high column densities and large covering 
factors. Although we cannot rigorously prove that intrinsic 


X-ray weakness is not present in these systems, the connec- 
tions of these quasars to other X-ray weak quasars includ- 
ing WLQs and super-Eddington accreting quasars point to 
a universal wind obscuration scenario for the weak X-ray 
emission found in type 1 quasars, or even type 1 AGNs in 
general. Multi-epoch X-ray observations of the intrinsically 
X-ray weak quasar candidates will further help clarify their 
nature. Besides variability investigations, deeper NuSTAR 
observations of PHL 1811 could provide further evidence of 
heavy X-ray obscuration. Also, higher signal-to-noise ratio 
and higher spectral resolution observations with future gener- 
ation X-ray observatories (e.g., Athena; Nandra et al. 2013) 
could reveal spectral features (e.g., the Fe lines in the repro- 
cessed component) that help discriminate between different 
scenarios. 
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APPENDIX 


A. XMM-NEWTON OBSERVATION OF LBQS 1442—0011 


LBQS 1442-0011 is a BAL quasar at z = 2.226 with a B-band magnitude of 18.2 (Gallagher et al. 2006). Its H6-based single- 
epoch virial SMBH mass is ~ 8 x 10°Mo, and its estimated Eddington ratio is zz 0.17 (Yuan & Wills 2003). Its previous Chandra 
observations have a co-added depth of 15.9 ks. Through systematic analyses of the Chandra observations of the Gallagher et al. 
(2006) LBQS BAL quasar sample, Liu et al. (2018) identified LBQS 1442 —0011 as one of the two good candidates for intrin- 
sically X-ray weak quasars based on its significant hard X-ray weakness (by a factor of fy = 1225) and its nominal hard X-ray 
spectral shape (Teff = 1.9*955, The other candidate is LBQS 1203+ 1530. 

Due to the large uncertainties of the fw and Teg values from the Chandra results, we proposed for deeper XMM-Newton 
observations of the two candidates, aiming to improve the parameter constraints. The targets were accepted at priority C, and 
LBQS 1442-0011 was observed on 2021 February 6 with a nominal exposure time of 87 ks. Unfortunately, the observation 
was affected significantly by background flares, and the cleaned exposure time is only 34 ks. Thus the sensitivity of the new 
XMM-Newton observation is only comparable to the previous co-added Chandra exposure. We processed the XMM-Newton 
data following the procedure described in Section 2.3. We chose a smaller source extraction region with a radius of 25" in order 
to increase the signal-to-noise ratio of this faint source. We also limited the upper-energy bound to 8 keV. The resulting net 
source counts are 48*15 in the 0.3—2 keV band (rest-frame 1.0—6.5 keV) and dd in the 2-8 keV band (rest-frame 6.5—26 keV). 
The T'effo.3-8 value inferred from the band ratio is 1.0487, We also fit the spectrum with a power-law model modified by Galactic 
absorption (phabs * zpow), and the best-fit I value (0.9 + 0.4) is consistent with the photometric result. The derived factor of 
hard X-ray weakness is fy 2 4 +2. 

Compared to the previous Chandra constraints, the XMM-Newton results suggest that the hard X-ray spectrum became flat- 
ter, and the observed hard X-ray emission became brighter (fw dropped). Therefore, we suggest that the X-ray weakness of 
LBQS 1442-0011 is also caused by variable partial-covering absorption, similar to the three quasars studied here (see Sec- 
tion 4.3 for discussion). 


